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Abstract 

In this paper, we study decoherence on Grover's quantum searching algorithm using a perturbative 
method. We assume that each two-state system (qubit) suffers cr z error with probability p (0 < p < I) 
independently at every step in the algorithm. Considering an n-qubit density operator to which 
Grover's operation is applied M times, we expand it in powers of 2Mnp and derive its matrix 
element order by order under the n — + oo limit. (In this large n limit, we assume p is small enough, 
so that 2Mnp(> 0) can take any real positive value or 0.) This approach gives us an interpretation 
about creation of new modes caused by a z error and an asymptotic form of an arbitrary order 
correction. Calculating the matrix element up to the fifth order term numerically, we investigate a 
region of 2M np (perturbative parameter) where the algorithm finds the correct item with a threshold 
of probability P t h or more. It satisfies 2Mnp < (8/5)(l — -Pth) around 2Mnp ~ and Pth — lj an d 
this linear relation is applied to a wide range of P t ^ approximately. This observation is similar to a 
result obtained by E. Bernstein and U. Vazirani concerning accuracy of quantum gates for general 
algorithms. We cannot investigate a quantum to classical phase transition of the algorithm, because 
it is outside the reliable domain of our perturbation theory. 



1 Introduction 

Since the idea of quantum computation appeared 000, a lot of researchers have been investigating its 
properties, algorithms, and implementations A quantum computer can be thought a sequence of 

operations which are unitary transformations and measurements applied to two-state systems (qubits). 
(The qubit means a system defined on a 2-dimensional Hilbert space {|0), |1)}.) For realizing perfor- 
mances that conventional (classical) computer hardly shows, it makes use of the properties of quantum 
mechanics, such as principle of superposition and its interference, principle of uncertainty, and entangle- 
ment (quantum correlation which is stronger than classical one) . 

One of the most serious problem for realizing quantum computation is decoherence, which is caused 
by an interaction between the system of quantum computer and an environment that surrounds it || H . 
It is pointed out that quantum information stored as a quantum state is fragile and collapses at ease 
by this disturbance. To investigate it, some decoherence processes are assumed and their effects on 
quantum algorithms are estimated [pl||. For overcoming these troubles, quantum error-correcting codes 



are proposed and their availability is examined 10 |ll| 



Not only for practical purposes but also for theoretical interests, it is an important question how 
robust the quantum algorithm is against this disturbance. We can expect that the quantum computer 
loses its efficiency gradually as decoherence gets stronger. Some researchers regard it as a quantum to 
classical phase transition Q . 

Grover's algorithm is considered to be an efficient amplitude amplification process for quantum states, 
so that it is often called a searching algorithm (l3) jl4j . By applying the same unitary transformation to 
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the state in iteration and amplifying an amplitude of one basis vector that we want gradually, Grover's 
algorithm picks up it from a uniform superposition of 2™ basis vectors with certain probability by 0(2™/ 2 ) 
steps. Because it handles a general problem (an unsorted database search), it can be formulated as an 
oracle problem, and it is proved that its efficiency is optimal in view of computational time (the number 
of queries for the oracle) |Q jL5] , many researchers have analysed this algorithm precisely and proposed 
a lot of applications jl6| . 

In this paper, we study the decoherence on Grover's quantum algorithm with a perturbative method. 
We assume a simple model and investigate it for higher order perturbation (numerically up to the fifth 
order correction), under the limit of an infinite number of qubits. 

The model has the following three characteristics. First, in the Grover's algorithm, we assume that 
we search the basis vector of 1 ■ - - 0} from the uniform superposition of the n-qubit logical basis. This 
assumption simplifies the iterated transformation. Second, we assume each qubit interacts with the 
environment independently and suffers a phase damping which causes tr z error with probability p and 
does nothing with probability (1 — p) jl7). Third, we take the limit of n — > oo, so that the matrix element 
of the density operator is simplified. 

In our perturbation theory, we expand an n-qubit density operator to which Grover's operation is 
applied M times in powers of 2Mnp . Investigating higher order terms of the perturbation, we obtain 
a physical interpretation that the a z error creates new modes as the algorithm goes steps. When we take 
the large n limit mentioned above, we assume p is small enough, so that a perturbative parameter 2Mnp 
can take any positive value or 0. Taking the n — > oo limit simplifies the matrix element of the density 
operator and gives us an asymptotic form of an arbitrary order term. 

Calculating the matrix element up to the fifth order term numerically, we investigate a range of 
2M np where the algorithm finds the correct item with a threshold of probability P t h or more. It satisfies 
2Mnp < (8/5)(l — Pth) around 2Mnp ~ and P t h — 1, and this linear relation can be approximately 
applied to a wide range of Pth- Hence, if we fix P t h to a certain value (Pth = 1/2, for example), we have 
to suppress the error rate to a value which is in proportion to the inverse of the number of quantum 
gates. 

Similar results are obtained by E. Bernstein and U. Vazirani in the study of accuracy for quantum 
gates They consider a quantum circuit where each quantum gate has a constant error because of 
inaccuracy, so it is an error of a unitary transformation and it never causes dissipation to the quantum 
computer. They estimate inaccuracy e for which the quantum algorithm is available against the fixed 
number of time steps T, and obtain 2Te < 1 — P t h- If we regard p/2 as inaccuracy e, and 2Mn as the 
number of whole steps in algorithm T, it is similar to our observations except for a factor. 

A. Barenco et al. study the approximate quantum Fourier transformation (AQFT) and its decoherence 
P|. Although motivation is slightly different from E. Bernstein and U. Vazirani's, we can think their 
model to be the quantum Fourier transformation (QFT) with inaccurate phase gates. They confirmed 
that AQFT can make a performance that is not so worse than QFT's one. 

This article is arranged as follows. In Section ||, we describe the model that we analyse in this paper. 
In Section [| we formulate a perturbation theory for our model and explain physical quantities that we 
estimate. In Section || we derive the matrix element of the density operator of the quantum computer for 
the 0-th and first order. We give a physical interpretation about creation of new modes by a z errors in 
Section ^|. Then, we derive the second order correction of the matrix element in Section |^. In Section ^, 
we take the limit of an infinite number of qubits and give the asymptotic form of an arbitrary order term. 
In Section @, we carry out numerical calculations of physical quantities up to the fifth order correction. 
In Section M, we give brief discussions concerned with our results. We collect formulas for deriving matrix 
elements in Appendix [a], and give some notes about numerical calculations of higher order perturbative 
terms in Appendix M. 



2 Model of decoherence 

In this section, we describe a model that we analyse. It is a quantum process of Grover's algorithm which 
suffers a phase error in iteration. 

At first, we give a brief review of Grover's algorithm Ga]. Starting from the n-qubit (n > 2) uniform 
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superposition on a logical basis, 



v ie{o,i}" 



it increases gradually an amplitude of a certain basis vector \xq) (xq £ {0, 1}") which is indicated by a 
quantum oracle. An operator W in Eq. (^) is an n-fold product of a one-qubit unitary transformation 
(Hadamard transformation) and given by W = H® n , where 

<0| <1| 

ff = tft i|o>(i M H t H = L (2) 



V2|i> V 1 - 1 

The quantum oracle can be regarded as a black box, and actually it is a quantum gate which shifts phases 
of logical basis vectors as 

R . / No) -> -No) /o>> 

rtxo ■ \ \x) -> |k) for z ^ x ' W 

where x ,x G {0, 1}", R Xo = R Xo , and Rl R Xo = I 

To let probability of observing \xq) be greater than a certain value (1/2, for example), we repeat the 
following procedure 0(\/2™) times. 

1. Apply i? Xo to the n-qubit state. 

2. Apply £> = WR W to the n-qubit state. 

Rq is a selective phase shift operator which multiplies a factor (—1) to |0 • • • 0) and does nothing to the 
other basis vectors, as defined in Eq. (|^). D is called the inversion about average operation. 

From now on, we assume that we amplify an amplitude of |0 - • - 0). From this assumption, we can 
write an operation iterated in the algorithm as 

DRq = {WR W)R . (4) 

After repeating this operation M times from the initial state of W\0) (= W\0 • • • 0) ), we obtain the state of 
(WR a ) 2M W\0). (We often write |0) as an abbreviation of the n-qubit state |0 • • • 0) for a simple notation.) 
Next, we think about the decoherence. In this paper, we consider the following one-qubit phase error 

p — > p' = p<r z pcj z + (1 - p)p for < p < 1, (5) 
where p is an arbitrary one-qubit density operator and a z is one of the Pauli matrices given by 

<0| <1| 
T t - 1°) ( 1 



= Ii^o _i| (6) 

For simplicity, we assume that the phase error of Eq. (^) occurs in each qubit of the state independently 
before every Rq operation during the algorithm. It assumes that each qubit interacts with its environment 
and suffers the phase error independently. 

Here, we add some notes. First, because Rq e U(2 n ) is applied to all n qubits and H € U(2) is applied 
to only one qubit, we can imagine that the realization of Rq is more difficult than that of W = H® n . 
Hence, we assume that the phase error occurs only before Rq. Second, although we assume a very simple 
error defined in Eq. (|^), we can think other complicated errors. For example, we can consider a phase 
error caused by an interaction between the environment and two qubits, and it may occur with probability 
of 0(p 2 ). In this paper, we do not assume such complicated errors. We discuss how the disturbance of 



Eq. (Bj) occurs in Section 9.2 
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3 The perturbative method 

Let /c/ M ) be the density matrix which is obtained by applying Grover's operation M times to the n-qubit 
initial state W|0). The decoherence of Eq. (||) occurs 2Mn times in p^ M \ We can expand p^ M ^ in powers 
of p and (1 — p) as follows, 

pi M ) = (1 - P f M ^ M) + (1 - pf Mn - x pT[ M) + ■■■ 

2Mn 

= Y,^-v? Mn - h M M \ (7) 



/i=0 



where {T^ M) } are given by 



T (M) = (WR ) 2M W\0)(0\W(R W) 2M , (8) 



2M-1 



T (M) = (WR ) 2M - k a^{WRo) k W\0)(0\W(RoW) k a^(RoW) 2M - k , (9) 

i=l fe=0 
n n 2M-1 

T 2 M) = EE E (Wiib^-VWaWO^JIoJ^lOXh.cl 



i=l 3=1 fc = 



n n 2M-12M-H 

EE E E (WR ) 2M - k - l a^(WR ) l 4 j HWR o ) k W\Q)(h.c.\, (10) 



and so on, cr£ represents the operator applied to the z-th qubit (1 < i < n), and (h.c.| represents a 
hcrmitian conjugation of the ket vector on its left side. We can reg ardli M) as a density operator whose 
trace is not normalised. It represents states where h errors occur during the iteration of M operations. 
On the other hand, we expand p^ M ^ in powers of p as follows, 

pM = p[ M) + 2Mnpp[ M) + l - (2Mn) (2Mn - l)p 2 p[ M) +■■■ 

2Mn , \ 

= e(TW m) > (") 



h=0 



where 



(M) _ T (M) 



AM) _ T (M) _,_ ?V 



(M) 



2Mn 

P1 M) = (-D^B-D^^y'Q)^ farfc = 0,l,...,2Mn. (12) 

In Section ^, we take a large n limit (the limit of an infinite number of qubits). If 2Mnp is small enough, 
we can consider the series of Eq. ( |ll| ) to be a perturbative expansion. Because we divide T^ M ' by (Mn) h 
as Eq. (jlj), an expectation value of pi (/i = 0, 1, • • • , 2Mn) can converge on a finite value in the limit 
of n — » oo. 

With these preparations, we will investigate the following physical quantities. Let Pth be a threshold 
of probability (0 < Pth < 1), so that if the quantum computer finds the item that we want (in our model, 
it is |0)) with the probability P t h or more, we regard it available, and otherwise we do not consider it 
available. Then, we consider the least number of the operations that we need to repeat for amplifying 
the probability of observing |0) to P t h or more for given p. We can describe it as M = M t h(p; Pth), and 
it satisfies (0|p( M )|0) = P t h- (For convenience, we write it as M th (p) with omitting P t h as far as it does 
not make confusion.) Because of 

(0|T (M) |0) P=0 - sin 2 (2M/V2") for large but finite n, (13) 
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M th (p)/2 



n/2 



(p=(l/2)arcsin (P t h) 




Figure 1: Variation of M t h(p)/v2" against p with the threshold probability P t h under large but finite n. 
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Figure 2: Variation of p c against Pth- 



obtained in Section [|, Mth(p) takes a value of <^v2™ at p = (with no decoherence) for large finite n, 
where ip = (1/2) arcsin VPh ■ 

As p gets larger, we can expect that M t h(p) increases monotonously. It could be possible for certain 
p c or more that we never observe |0) at least with probability of P t h- (Hence, p c depends on Pth-) Such 
behaviour of M t h(p) can be drawn in Figure [l]. We multiply a factor to M t h{p) for normalisation. 

Because (0|p( M ) |0) with p = increases monotonously from M = to M = (7r/4)v2" and then decreases 
as shown in Eq. (|l3|), M t h(p)/v / 2" at p = p c is equal to or less than 7r/4. 

Regarding P t h as a threshold whether the quantum computer is available or not, we can consider p c 
to be a critical point. (This is not a so-called quantum to classical phase transition.) We can draw a 
graph of p c against P t h as Figure |[ (If P t h = 1, we obtain p c — 0.) 

In this paper, we calculate physical quantities using the perturbative parameter x = 2Mnp, so that 
we take M/\/2™ and x for independent variables. (In our original model defined in Section |^, we take 
M and p for independent variables.) We can define as well M t h(^; Pth)/\/2™ that represents the least 
number of the operations iterated for amplifying the probability of |0) to Pth for given x. Furthermore, 
we also obtain x c , for which or more we can never detect |0) at least with probability P t h- Hence, we 
obtain a graph of Mth(^) / V^2" versus x instead of Figure 0, and that of x c versus Pth instead of Figure ||. 
The differences of these quantities are discussed in Section [|. 

The dependence of x c on P t h gives us useful information. If we regard 2AIn as the number of 
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computational steps T, and p/2 as a parameter e which represents a degree of errors, it serves the region 
of Te where the quantum computer is available for P t h, because of x — 2Mnp. To make these analyses, 
we need to know (0|p^ M -'|0). In the following sections, we calculate (0|/3^ M ^|0) to the fifth order of p (up 
to (2Mnp) 5 ) numerically under the n — > oo limit. 



4 Matrix elements up to the first order 



7T = [sin(2M + 1)*|0) + COS( ; 9 :_ 4 - 1)g X>)](h.c.|, (14) 



In this section, we consider matrix elements of the density operators with no and one error, (0|Tq M ^|0) 
and (0\T^ M) \0), defined in Eqs. (§) and (|). 

First, we derive T (M) and (0|T (M) |0). Using Eq. @ in Appendix [vG], we obtain 

<M)_^m , ^ fl , n y , cos(2M+l)0 

(0|T (M) |0) = sin 2 (2Af + 1)0, (15) 

where 

1 /2™ — 1 
sin0 = ■ , cos0 = \ . (16) 

V 2™ 

(This parameter is introduced by M. Boyer et al. and it simplifies our notation [H.) From Eq. (p"5j), 
we notice the following facts. If there is no decoherence (p = 0), we can amplify probability of observing 
|0) to unity. Taking large (but finite) n, we obtain sin0 ~ and ~ 1/V2™, and we can observe |0) with 
probability of unity after repeating Grover's operation Af max ~ (7r/4)\/2™ times. 

Then, we think about t[ M ^ and (0|t| m ^|0). For convenience, in spite of the definition of t[ M ^ given 
in Eq. ([)]), we rewrite it as follows, 

n M-l 



T i = ^(^o) 2(M " fe V«(IFi? ) 2fc W / |0>(h.c. 



i=l fc=0 
ri M-l 

+ Z Z (^o) 2(M " feK Vi 4) (IUi 1 , o) 2 ' £+1 VF|0)(h.c.|. (17) 

z=l fc=0 

We can derive an explicit form of (0|Tj |0) using formulas collected in Appendixes Al and A. 2. The 
matrix element of the first term in Eq. ( |l7j ) can be given by |7^icn| 2 , where 

T e W = {0\{WR f( M -Vaf{WR Y k W\Q) 



= (-!)-"* [cos 2(M- k)9(0\ + Sin2 /if~ fc)g E^0^ ) 



x(-l) fc [sin(2fc + 1)*|0> + ^J=^E I»>1 

(-1) M [cos 2(M - fc)0 sin(2A: + 1)0 — sin 2(M - fc)0 cos(2fc + 1)9], (18) 



1 



We notice that 7~Jmn does not depend on the subscript i. In a similar way, we can obtain a matrix element 

(*0 12 



of the second term in Eq. ( |l7| ) as |^ d( j| 2 , where 



= (-l) M [sin(2M - 2/c - l)0cos2(fc + 1)0 — cos(2M - 2k - 1)9 sin 2{k + 1)9}. (19) 

2" — 1 

(Here, we notice |7^ V cn| 2 = \T^~ k ~ 1 ^ | 2 -) Hence, we obtain 

M-l 

(0|lf>) = n E d^nl 2 + (l^dil 2 )- (20) 

k=0 



5 Physical interpretation for creation of modes 



If we derive the explicit form of T^ M , we can obtain a physical interpretation of creating modes by a z 
error. 



Let us consider the state included in Eq. (|17[), 

(WR o ) l a^{WR o ) k W\0}, 



(21) 



which suffers the phase error only once. From Eq. (|99|) in Appendix A.l, we obtain the following equation, 
a ( i ) {WR ) 2k W\0) 

= (-l)*[ S m(2fc + 1)0|O) + C °^ + | )0 (E l*> ~ E I-))] 



V2^ 

(WR ) 2k W\0) - V2(-i) k COs{2k+ n 1)e \ n 



7i» 
= 



for k = 0, 1, 



cosf 



where 



V2 



Eire) for i = 1, • • • , i 



In a similar way, from Eq. (IOC), we obtain 



a^(WRo) 2k+1 W\0) = (WR ) 2k+1 W\Q) + v^(-l)' 



,sin2(fc+ 1)0, 
cos 6* 



I* 



for jfc = 0, 1, 



(22) 
(23) 

(24) 



To obtain an explicit form of Eq. (|2l|), we have to apply WRq to Eqs. ( |22j ) and ( |24| ) from the left side 
step by step. (We count WRq as one step for a while.) We can derive an explicit form of (WRo) l \r]i) as 
follows. In the case of I = 1, we obtain 



WR \r,i) = -^=(|0) - |S)), 

where i represents a binary string whose all digits are '0' but the i-th digit is '1', so that 

\i) = jO — 1 0---0). 

T 



(25) 



(26) 



For I = 2, we obtain 



Therefore, we obtain 

(WR ) 2l K 
(WRo) 21 ^^ 



(WR a ) 2 \ m ) = ±WR (\0)-\i)) 

= Ei*> 

a:i=0 
V^WjO) + 1 77;}. 



i-1 



-V2 ^ (VFi? o ) 2m T^|0) + \rn) for Z = 0, 1, • • •, 

771=0 
l-l 

-V2 (WR o ) 2m+1 W\0) + -=(|0) - |i» for I = 0, 1, • • 

m=0 * 



(27) 

(28) 
(29) 



where Y2i=o means that no term is summed up. 
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W|0> - 
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II 



2(k+l) 
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Figure 3: Creation of modes by the error. 

From Eqs. (§|), @, ©, and ©, we can describe (W-Ro)'^ (W.Ro)*W|0) for fc = 0, 1, • • • and 
I = 1,2, ••• with (VKi? o ) fe W|0), and (|0) - (We obtain the explicit form of (WR o ) k W\0) in 



Appendix Al.) For example, we can write 
(WR ) 2l <r[ l) (WRo) 2k W\0) 

= (WRo) 2(k+l ^W\0) - V2(-l) k COS{2k + 1)0 [-V2 £ (WR ) 2m W\0) + fa)] 



m=0 

for k = 0,1,- ■■ and I = 1,2,- • -. (30) 

This equation allows us the following interpretation. If a z error occurs in the i-th qubit of the n-qubit 
state at the 2fc-th step, it causes new modes which are created as the initial state W\0) at every two steps 
from 2(k + 1), that is, 2(k + l)-th, 2(k + 2)-th, 2(k + 3)-th, • • •, and so on. (See Figure |.) The state 
becomes a superposition of them. 

Here, we derive the matrix element (0|t| M ^|0) again using this interpretation. From Eqs. ( |30"| ) and 



in Appendix A.l, we obtain 



(0\(WRo) 2(M - k) a^(WRo) 2k W\0) 
= (0\(WR ) 2M W\0) 

-72(-l) fc C ° s(2fc ± 1)6 [-V2 M ^ (0\{WRo)*W\0) + (0\ Vi )\ 
cos — ' 

l=o 

= (~1) M sin(2M + 1)9 + 2 (-l) fc COs(2fc + 1)0 \-l)' sin (2Z + 1)0. (31) 

cosfl ^ 



(We eliminate one a z operator from the above equation.) Then using a formula of Eq. (103) in Ap 



pendix A. 3 to sum up trigonometric functions, we can rewrite Eq. ( |3l| ) as 
(0| {WR Q ) 2{M - k) ai l) (Wi2o) 2 *W|0) 

= (-1)" sin(2Af + 1)0 + 2(-l)^ + 1)6 ( - 1)M '" 1 ^ 2(M = k)d 
= (-l) M [cos2(M - fc)(9sin(2fc + 1)6* - — - - sin2(M - k)6 cos(2fc + 1)9} 

where we substitute cos 2 6 — (2™ — l)/2™ of Eq. (pf). In this derivation, we expand the matrix elements 



into a series of modes by Eq. (|30|) , and sum up them by the formula of Appendix [A.3| . We often use this 
technique in this paper. 
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Figure 4: Relation between S, S%, and <S 2 - 



In a similar way, using Eqs. (Um, (|2^), ( |100| ), and (10£), we obtain 

(OliWRof^-^cr^iWRof^WlO) 
= (0\(WR a ) 2M W\0) 



M-k-2 



^ ( _ i)fc sin 2^+1)^ £ (0|(Wflb) 9,+1 W|0) + -L(0|(|0>-|i»] 



= (-1)* sin(2Af + 1)0 + 2 (-l) fc Sin2(fc + 1) V £ '(-I)' cos2(/ + 1)9 + h 

COS u 

1=0 

= (-1) M sin(2M + 1)9 

+2 (-1 ) k v n ' [-- - ±— '- — cos (2M -2k- 1)9 + - 

v ; cos6» 1 2 2cos6< v ; 2 J 

= ^.dd- (33) 



6 Matrix element of the second order 



In this section, we consider the matrix element (OlT^^jO) which contains two a z errors, defined in 
Eq. (|l0|). We make good use of the interpretation of creating new modes discussed in Section || for 
obtaining it. 

Let us see the first term of Eq. (lOh. It suffers two a z errors at the same step as follows, 



(0\(WRo) 2M - k a^4 j) (WRo) k W\0) 



for 1 < i < n, 1 < j < n, and i =/= j . 



(34) 



Here, we consider the following term, 



{Q\{WR Q )^ M -^afa{^{WR Q ) 2k W\Q) 
= (-l) M [cos2(Af - fc)6>sin(2fc + 1) 



1 



2™ - 1 
for k 



sin2(M- k)9 cos(2k 
0,l,---,Af-l, 



x^Oy^O 



(35) 



where we use Eqs. ( |99] ) and (101) in Appendixes A.l and [A.2| . To obtain an explicit form of the above, 

we have to calculate J2 x ^o ^2 V jto( x \ a ^ a ^ \v) for i ^ 3- 

For deriving it, we define a set S = {\x) : x £ {0,1}", x ^ 0}, and its subsets, S\ = {\x) : \x) £ 
S, Xi — 1} and 52 = {\x) : \x) £ <S, Xj — 1}, as shown in Figure |]. The number of elements of them and 
S\ n S2 are given as |<S| = 2™ — 1, |5i| = ]«Sa | = 2" _1 , and |<Si n S%\ = 2 n ~ 2 . A set of basis vectors whose 
signs are flipped by ai is given by (6>i U S2) H (Si PI S2), and the number of its elements is equal to 
2 ra_1 . Hence, the number of elements in S whose signs are not flipped is equal to (2™ _1 — 1). From these 
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considerations, we obtain 



and 



(0\{WR ) 2 ( M -V<T^ai j HWRo) 2k W\0) = T e « 
In a similar way, we can obtain 



(OKWRo^-^a^a^iWRof+'WlO) = T o {k) 



del 



for k = 0, M-l. 



(36) 
(37) 
(38) 



Therefore, the matrix element of the first term in Eq. ( |l0|) is equal to [(n — 1)/2](0|T^ M ^|0). 

Then, we think about the second term in Eq. (|l0|). Using Eq. (pG), we obtain the following term, 



rj-(k,l;Sij ) 
even, even 



(0\(WRo) 



2( M -k-l)(i) 



{WR ) 2l a^(WR ) 2k W\0} 



(0\{WR Q ) 2{M - k - l) a^(WR ) 2{k+l) W\0) 
_^2 ( _ 1)fe cos(2fc + l)0 



cos 8 



i-i 



x[-V2J2 {0KWRo) 2{M - k - l) a^(WR o ) 2m W\0) + {0\(WR a ) 2 ( M - k ~ l) a^\ Vj )}. (39) 

The bracket [• • •] in the first line of Eq. ( |39| ) represents that this part is calculated at first. Here, we use 
the same technique in Section [5] again. We expand the matrix element by modes caused by the a z error 
and sum up them. (We eliminate one a z operator from Tevdn^even-) As a result of Eq. (|39j), we obtain 
terms which contains only one a z error, and essentially they have been obtained already in Sections ^| 
and |^. Here, we introduce a notation of 



gW(k,l) = (0\{WR ) l a^(WR ) k W\Q), for k = 0, I = 1,2,- 



(40) 



and collect its explicit form in Appendix A. 4. We also collect some formulas of \r)j) and (|0) — \j}) in 
Appendix A. 5. 



Using Eqs. ( |103| ), ( |105| ), ( |107| ), and ( pLllD in Appendixes |AJ3|, |AJ, and|Aj>| we obtain 

7±^l = G^(2(k + l),2(M-k-l)) 

k cos(2k + l)8 



even, even 



-V2(-1) A 



cos 8 



: {-V2^(-r 



M-k-l+m 



[cos 2(M -k-l)6 sin(2m + 1) 



m=0 



1 



sin 2(M -k-l)6 cos(2m + 1) 



2" - 1 

(-l)M-k-i s [ n 2(M - k - 1)8 



V2 cos6» 
= g^(2(k + l),2(M-k-l)) 



-(-1) 



M cos(2/c + 1) 



cos 2 8 



x { [cos 2(M -k-l)6 sin 2Z<9 - — - sin 2(M - k - 1)8 cos 2Z6>] 

+(- 1 )" i (^T3y - <%) sin2(Af - fc - 1)8}. 
From the above, we notice that 7^ion',even depends on (and not on i and j). 



(41) 
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As results of similar considerations, using Eqs. (|22|), (E4), (|28|), and formulas in Appendixes A. 3 



A. 4, and A. 5, we obtain the other matrix elements 

(k,l;6 i} ) 



odd, even 



r 



even, odd 



^j-(k,l;Sij) 
i odd,odd 



(0|(T^i?o) 2(M_fe_0_1 ^(VKi? ) 2 ' ( j( j) (V^i?o) 2fe+1 W"|0) 

g^(2(k + l) + l,2{M-k-l)-l) 
M sin2(fc + 1)0 



cos z 



x {[sin(2M -2k -21 -1)0 sin 219 + cos(2M - 2k - 21 - 1)9 cos 2Z0] 

2™ — 1 

+ - %) cos(2M -2k -21- 1)9}, 



(42) 



(k,l;Sij) _ 



= (0\(WR ) 2 < M - k - l >- 1 v® 



(WR ) 2l+1 ai f) {WR ) 2k W\0) 



G^(2{k + l) + l,2{M-k-l)-l) 



-V2(-l) 



k cos{2k + l)9 
cos 9 



i-i 



-V2 G {1) (2m+1,2(M -k-l)-l) 



+ _L ( o|(^i?o) 2(M - fc - () -VW(|o)) - \j)}} 

G {1) {2{k + l) + l,2(M-k-i)-l) 
-V2{-l) kCOa{2k + 1)9 



cos 9 



:{-V2£(-l) 



M-k-l-l+m 



[- sin(2M -2k -21- 1)9 cos 2(m + 1)0 



m— 



•— — - cos(2M -2k -21- 1)9 sin 2(m + 1)0] 



(-1) 



M-fc-Z-l 



-[sin(2M -2k -21- 1)9 + (-1)' 



V2 

S (1) (2(A; + + l,2(M-fc-0-l) 
M cos(2fc + 1)0 
1 j cos 2 

x{[sin(2M -2k -21- l)9cos(2l + 1)9 - 



. cos(2M -2k -21- 1)0 
V2" - 1 



]} 



■ cos(2Af -2k -21 -1)9 Bin(2Z + 1)0] 



(-l)'r 1 



(-!)*«] cos(2M -2k -21- 1)0}, 



V2" 2™ - 1 

(0\(WR )^ M - k - l -^a^\WRo) 2l+1 a^\WR ) 2k+1 W\0) 

g^(2(k + l) + l,2(M-k-l-l)) 
M sin2(fe + 1)0 



(43) 



cos z 



x{-[cos2(Af -k-l- l)9cos{2l + 1)9 + 



2™ - 1 



sin 2(M -k-l -1)9 sin(2Z + 1)9} 



(-l)'r 1 



^2™ 2™ - 1 
write (0|T 

(0|T 2 (M) |0) 



+ (-l) 5 ^] sin2(M -k-l- 1)9}. 



(44) 



Finally, we can write (OlX^^lO) as 



n-1 



(OIT^IO) 



M-l M-k-l 



^ t ^ , [fl(fi 1) l^cven',evenl + n l^cvc'n',evenl 



fc=0 1=1 
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+ E E Kn-l)|^ n | 2 + n|^ n | 2 ] 
fe=0 1=1 
M-l M-fc-1 

+ E E Kn-^lT^^P+H^Ll 2 ] 
fc=o ;=o 

M-l M-fe-2 

+ E E K»-i)i^S2di a +»i^S!2di a ]- (45) 

fe=0 /=0 



7 Large n limit and asymptotic forms of matrix elements 

The matrix elements, (0\T^ M) \0) and (0|T 2 (M) |0) obtained in Sections | and ^, are too complicated to 
handle as they are. In this section, we take the limit of an infinite number of qubits (n — > oo), and 
discuss their asymptotic forms. We also discuss how to obtain an asymptotic form of any higher order 
term under n — ► oo. 

We consider the limit of n — > oo for n-qubit state. We assume we can take very small p, so that 
x = 2Mnp can be an arbitrary real positive value or 0. If < M < 2n^/2^, M0 converges on a certain 
value of 8 (0 < 8 < 2ir) under this limit. (The definition of 9 is given in Eq. fllfij).) It is reasonable that 
we assume M is order of 0( V / 2") or less and define 6 = limbec M9. Because (0\T^ M) \0)/(Mn) k does 
not depend on p, we can take n — > oo for it naively. Hence, from Eq. (|l5|), we obtain 

lim (0|T (M) |0) = sin 2 28. (46) 

n — >oc 

Then, we consider the asymptotic form of (0|T^ M ^|0) in Eqs . (|l8| ), (pj|), and (20). To let it converge 
on finite value, we divide it by a factor of Mn as shown in Eq. (|12|). We can obtain 

(Qljf^lQ) u i_ y^nf(fc) i2 + ,^(fe) ,2i (47) 

fc=0 

where 

= (-l) M cos2(Af -fc)0sin(2A: + l)e, 
f ( d ^ = (-l) M sin(2M-2fc- 1)6> cos 2(A: + 1)9. (48) 

(We drop the terms with a factor 1/(2" — 1) in and and obtain Eq.(E§).) We substitute 

8 = limn-^oo M0, cj> — k9, and 

M-i 

lim V 9 = / d0 (49) 

j n U 



fc=0 



into Eq. fl47|), and obtain 

llf^O) 



lim A^U W = / d0{[cos2(e-0)sin20] 2 + [sin2(e 

n— ioo Mn 8 Jq 



cp) cos2</>] 2 } 



i _ 1^840-^^46. (50) 



Next, we consider an asymptotic form of (0|T 2 (M) |0) obtained in Eqs. @, (§|), @, @, and @. 
Because of convergence, we divide it by (Mn) 2 as Eq. ([l2|). In the limit of n — > oo, we can neglect 
(OIT^IO) and T Q ( ^ ;1) for a,0G {even, odd}, and we obtain 

/nl^( M )ln\ i M-l M-fc-1 M-XM-k-X 

i irn A^±2 m _ lim _L[V V if( fc ^°) l 2 +V V i-f (fc ' Z;0) I 2 

v ; fc=0 1=1 fc=0 1=1 

M-XM-k-X M-XM-k-2 

+ E E l^evc'n.oddl 2 + E E l^odd,odd| 2 ]' (^1) 
fe=0 1=0 fe=0 Z=0 
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where 



qr(k,l;0) 



even, even 



r 



(fe,i;0) 



odd, even 



r 



(k,l;0) 



even, odd 



V odd.odd 



(-l) M [cos2(Af - k - l)9sm(2k + 21 + 1)8 

- cos(2fc + 1)8 cos 2{M - k - 1)9 sin 2Z6>] 
-l) M cos2(M- fc- Z)0 cos 2Z0sin(2fc + 1)8, 
-1) M sin(2M - 2fc - 21 - 1)8 cos 218 cos 2(fc + 1)0, 

-l) M [sin(2Af -2k -21- l)0cos2(fc + Z + 1)0 
- cos(2fc + 1)6» sin(2M -2k -21- 1)8 cos(2Z + 1)0] 
-1) M+1 sin(2M -2k -21- l)0sin(2Z + l)0sin(2fc + 1) 
- I - 1)0 sin(2Z + 1)0 cos 2(fc + 1)0. 



l) M cos2(M 



(52) 



(We use lim jwoo cos0 = 1 because of Eq. (|l6|).) This asymptotic form contains only the terms where 
a z errors occur at different steps and at different qubits from each other. Hence, defining ip — 19 and 



lim^oo J2ilo k = Jo 9 dip, we obtain 



e-0 



lim 



(Q|^ (M) |Q) 

{Mn) 2 



1 

Q2 



d<p 



e-4> 



dip 



{[cos2(8 — <ft — ip) cos 2ip sin 20] 2 
+ [sin 2 (8 - - ip) cos 2y> cos 20] 2 
+ [sin 2(8 — 4> — ip) sin 2<p sin 20] 2 
+ [cos2(8 — 4> — ip) sin2tpcos20] 2 } 
3 



i-^cos48 



648 



sin 48. 



(53) 



Seeing Eqs. (48), (f50|), (p2|), (|53|), and formulas of Appendix A. 4 , we find how to obt ain the as ymp totic 
form of h-th density operator (h = 1, 2, • ■ •) under n — > oo. We derive it in Appendixes |A.6|, A. 7, and A. 8 



Here, we use only its result. Preparing an /i-digit binary string a — (cti, ■ ■ ■ , ah) G {0, l} h , we define the 
following 2 h terms, 



(201 



sin 
cos 



for h 



cos 
sin 



(202 



COS 

sin 



(20/,) 



cos 
sill 



(2(8 



h 

£< 

s=l 



1,2,. 



(54) 



where 



(,:) 



/(a;) for a = 
g(x) for a = 1 



We notice that the function of 0i and the other functions of 02, • • •, <fih, 8 — 52 s=1 
and cosine functions are put in reverse). These terms are integrated as 



lim 

n—>oo 



(Q|jf f) |Q) 

(Mn) h 
1 

Qh 





r e-4>i 








01 / 


02 •• 


L 


10 


Jo 







0-4,! 



h ^ \T au ..., ah [ 

( ai ,---,a h )<E{0,l} h 



(55) 

are different (sine 



',0h)| 2 . (56) 



We can find that the expression of Eqs. (|54|) and (|56[) is coincide with Eqs. (|50|) and (]5 



Here, let us calculate the asymptotic form of (0|T 3 |0). From the above rules, we obtain 



lim 

n — >oo 



(Q|jj M) |Q) 

(Mn) 3 



1 

Q3 





/•0-4>i 


/•& — <pl—4>2 






01 / a 


02 / a 


03 


10 


Jo 


Jo 





13 



{[sin 20i cos 202 cos 203 cos 2(9 — 
+ [cos 20i cos 202 cos 203 sin 2(0 - 
+ [sin20i sin 202 cos 203 sin 2(0 — 
+ [cos 20i sin 202 cos 20 3 cos 2(0 - 
+ [sin20i cos 202 sin 203 sin 2(0 — 
+ [cos 20i cos 202 sin 20 3 cos 2(9 - 
+ [sin 20i sin 202 sin 203 cos 2(0 — 
+ [cos 20i sin 202 sin 203 sin 2(0 — 

1 + 320 



1 
12 



ft 2 



3-160 
15360 2 



cos 40 



2O480 3 



01 — 02 — 


h)f 


- 01 — 02 — 




01 — 02 — 


h)f 


- 01 — 02 — 


h)] 2 


01 - 02 ~ 


h)\ 2 


- 01 - 02 - 


h)} 2 


01 - 02 ~ 


h)} 2 


01 - 02 - 


h)} 2 } 


2 




- sin 40. 





We pay attention to the following facts. If we expand the asymptotic forms of Eqs. (p 
and rt57n in powers of 0, we obtain 



lim 



3|^ M) |0) 
(Mn) h 



40 2 + 0(0 3 ) 
(8/3)0 2 + O(0 3 ) 
O 2 + O(0 3 ) 
(4/15)O 2 + O(0 3 ) 



for h 
for h 
for h 
for h 



(57) 



(58) 



Hence, they converge to under the limit of — > (or M — > 0). This means that the probability of 
observing |0) for the uniform superposition is almost 0, and it is reasonable. 



8 Numerical calculations of physical quantities 



In this section, we carry out numerical calculation of ( 



JM) 



|0) for the asymptotic form under the n — > oo 



limit, and investigate physical quantities explained in Section |3|. Especially, we discuss the critical point 
x Cl over which the quantum algorithm comes not to be available for the threshold probability Pth- 
In the perturbation theory, we can rewrite Eq. (fil|) under the limit of n — > oo as 



Prob(0,a;) 



lim 

n— ►oo 



->( M ) i 



pv~'|0) 

C o (O) + Ci(0)x- 



-C 2 (0)a; 2 



oo 1 



h=0 



hi 



(59) 



where O = lim„_ 



M0, x = 2Mnp, 

C (O) = F (Q), 
Ci(0) = -F (6) 



-Me), 



C h (Q) = (-^tjl(-\yj^-y F^) for h = 0,1, 



3=0 



(60) 



and 



F h (Q)= lim 



(0\T^\0) 



for h = 0,l,---. (61) 
(H), and @. Using the rules of Eqs. (|J) and 



(Mn) h 

FfJQ) for h = 0, 1, 2,3 are obtained in Eqs. ©, (fo 
( p6| ) , we obtain higher order terms in Appendix [b|. 

The original model that we define in Section || has two independent parameters, M and p 7 and n 
takes a fixed finite value. On the other hand, the representation of Eq. ([39]) has O and x as independent 
parameters, and n gets infinity, so that it does not have a certain fixed value for n. (It is clear that 
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0.5 







x 







0.5 



1 



1.35 



Figure 5: Variation of P ro b(0,:r) against x with fixed = 7r/4. (It means M is fixed to M, 



(7r/4)\/2™.) A thin dashed curve, a thin solid curve, and a thick solid curve show perturbations up to 

the first, third, and fifth order each. Black circles represent results obtained by Monte Carlo simulations 
of n = 8 case (8 qubits) with M max = 12. Each circle is obtained for x — 2M max np = 192p, where p is 
varied from 5 x 10 -4 to 5.5 x 10 -3 at interval of 5 x 10~ 4 . In these simulations, we make 20000 trials for 
taking an average. 

and x are independent of each other from their definitions.) These difference reflects physical quantities 
of M t h(p, Pth), Pc, M t h(x , Pth) , and x c , which are introduced in Section ^| When we estimate these 
quantities numerically, we examine their meanings and differences. 

In this section, we make numerical calculations up to the fifth order correction. To investigate the 
range of x where our perturbative approach is valid, we need to estimate the sixth order term of Eq. (|59|). 
From the discussion in Appendix [b|, we can consider it is reliable around < x < 1.35, so that the sixth 
order correction of P ro b(©,a;) is bounded to 10~ 3 . We compare the perturbation theory with results of 
Monte Carlo simulations of our model, and confirm its reliability in Figures]^ and || In these simulations, 
setting n = 8 (8 qubits), we fix p and cause <j z errors at random on each trial. We take an average of 
(0|p (M) |0) p , the probability of observing |0) at the M-th step (M = 0, 1, • • • , M max (= 12)), with 20000 
trials for each certain value of p. 

In Figures [7] and we consider perturbations up to an odd order (the first, third, and fifth). If we 
sum up even number of correction terms, (0|/9^ M '|0) with fixed (or M) does not decrease monotonously 
against x = 2Mnp (or p) . It turns for increasing from some value of x (or p) , so that sometimes we cannot 
find p c (or x c ) by numerical calculation. Hence, we always consider corrections up to an odd order. 
Figure |^ shows a variation of P ro b(0, x) against x with fixed O = 7r/4, namely M is fixed to M max = 



(it / 4)^/2~^ . (Hence, the independent parameter is only p actually, but n is infinite.) We draw the curves 

of Eq. d59| ) up to the first, third, and fifth order corrections, and plot simulation results. At x — 0, there 
is no error and P ro b is equal to unity. As the error rate x gets larger, P ro b decreases monotonously. 

Figure]! shows a variations of P ro b (©, x) against with fixed p. Because we use the variable x — 2Mnp 
instead of p in the perturbation theory, we have to rewrite 



and give some finite n. In Figure |6|, we set n = 8 and draw curves of perturbation theory up to the fifth 
order against with fixed p (p = 2.0 x 10~ 3 for the thick solid curve, and p = 5.5 x 10~ 3 for the thick 
dashed curve). We also plot results of the simulations. 

From Figure ^, we notice that the maximum value of P ro b is taken at O < 7r/4 for each p (we show 
these points with vertical thin dashed lines), and the shift gets larger as p increases. It means 8th (Pc, Ph) 
gets smaller than 7r/4, as P t h decreases. (We write O t h(p;Ph) = lim^oo M t h(p; Pth)#, and M t h(p;Pth) 






(62) 
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Figure 6: Variation of P ro b(0,x) (perturbation up to the fifth correction) against with fixed p. The 
thick solid curve and the thick dashed curve represent p = 2.0 x 10 -3 and p = 5.5 x 10 -3 each, for n = 8. 
Black circles are results of Monte Carlo simulations. These two curves show ©th(Pc) is lower than 7r/4. 
It is confirmed by the simulation results. 



means the least number of the operations iterated for amplifying the probability of |0) to Pth under the 
error rate p.) 

Then, let us see the behaviour of the algorithm with fixing the threshold of the probability on P t h- 
Figure represents a variation of 8th (p) against p with n = 8 (8 qubits) for P t h = 1/2. Seeing it, we can 
confirm that Oth (p) cannot reach to tt/4, even if p — p c . It is consistent with results of Figure^. 

When we draw curves of Figures g and 0, we have to put finite positive n. (We set n on 8.) This 
treatment cannot be fully justified, because Eq. ( |59| ) is obtained with the n — > oo limit. Next, we 
compute physical quantities with taking independent parameters © and x. (We need not give finite n.) 
Figure |] shows a variation of 6th(^) against x with P t h = 1/2, where ©th(£; Pth) = ^co-n-*oo M t h{x] Pth)#- 
(Mth(x; Pth) represents the least number of the operations to amplify the probability of |0) to Pth under 
given x.) Seeing Figure 0, we find &th(x) increases as x gets larger from x = 0, and it reaches to the 
maximum value at x = x c . Comparing Figures [7] and 0, we notice &th(x c ) > Oth (j>c) for P t h = 1/2. (We 
can actually confirm ©th(^c) > ©th(Pc) for < VPth < 1 by numerical calculations.) 

It can be explained as follows. Figure || illustrates a variation of (0\p( M ^\0) against O for < p < p c . 
In the case of p = (no decoherence) , we can obtain ©th(£ = 0) = ©th(p = 0) = ir/8 for P t h = 1/2. (If 
p = 0, we obtain x — 2Mnp = 0.) If we let p get larger until p = p c , ©th(p) increases gradually. Through 
this process, both x and t h(2;) increase. When we reach at p = p c , we obtain t h(a;) = ©th(p c ) < 7r /4 
for x = 2M th (p)np c < x c . 

Although p gets the allowed maximum value of p c , we want to increase both x and @th(£) still more. 
To increase x, we take the following trick. We decrease p by infinitesimal Ap, as shown in Figure ^|. 
Then, (0|p (M) |0) takes P t h = 1/2 at two points of O, and we write them as ©1 and 2 (©i < © 2 ). At 
this time, we take the large one of them as ©th(^), so that ©th(a;) = ©2. 



Because 



JM) 



|0) takes the local minimum value at ©th(Pc) for P = Pc, we can obtain 



d(0\ P ^\0) 



<90 



d(o|p( M )|o) 



e=e th ( Pc ) 

P = Pc 



d(rV M >|o) dp 



and 



dp 
06th (p) 



09th (p) 

c)(0|p( M )|0) 
dp 

= or 



P=Pc 



dp <9©th(p) 



P=Pc 



^0, 



dOth(p) 
dp 



(63) 
(64) 
(65) 
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Figure 7: Variation of 6th(p) against p with fixed n(= 8) for P t h = 1/2. A thin dashed curve, a thin 
solid curve, and a thick solid curve show perturbations up to the first, third, and fifth order each. The 
algorithm cannot observe |0) with the probability of 1/2 or more for p > 0.00589. Hence, p c w 0.00589 is 
the critical point of P t h = 1/2. We obtain 6 t h(p c ) ~ 0.2137T < 7r/4. 




Figure 8: Variation of Qth(x) against x with Pth = 1/2. A thin dashed line, a thin solid line, and a thick 
solid line show perturbations up to the first, third, and fifth order each. The algorithm cannot observe 
|0) with the probability of 1/2 or more for x > 1.12. Hence, x c w 1.12 is the critical point of P t h = 1/2. 
We obtain Gth^c) « 0.247tt. 
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cOlp (M) IO; 



P th =l/2 




©th(Pc) 



Figure 9: Variation of (0|/9 (M) |0) against 9 with < 
Pth = 1/2- White circles represent 0th (p) and &th(x). 



P < Pc- The threshold of probability is set to 



Hence, the difference of [02 — ©th(Pc)] can be quite large, and &2{p c — Ap) > &th(p c )Pc is possible. 
Therefore, we can make ®th(%) > ©th(Pc) and a; > 2M t \ l (p c )np c - (These considerations can be applied 
to < VPth < 1 as well.) 

Then, we move on to the variation of x c against P t h, which is shown in Figure [l0| We obtain it 
as follows. We calculate 0th(aO for given P t h as varying x from 0. (We use the Newton's method for 
obtaining a root of for the equation of P ro b(07a;) = Pth with given x.) When x gets a certain value, 
we cannot find a root for t h(a;), and we regard it as x c . By repeating this calculation, we obtain the 
curve of Figure [l^. 

In these calculations, we notice 0th(^c! Pth) — 7r/4 for the range of x and Pth where the perturbation 
theory is reliable (0 < x < 1.35). It is caused by the approximately symmetric property of P^(0) obtained 
in Section [?| and Appendix [b| as 



F h ((n/4) + A) ~ F h ((ir/4) - A) for < A < (vr/4) and h = 0, 1, 

However, strictly speaking, 0th(^cjPth) cannot be a constant for Vxc and VPth. 
Using Eq. ([59|), a tangent at Pth = 1 is given by 



c(l - Ph), 



1 



Ci(tt/4) 



(66) 



(67) 



because ©th(^c) — tt/A and x c — for P t h = 1- It means that the algorithm is available for 2Mnp < 
(8/5)(l — Pth) around Pth — 1, and this relation approximately holds for a wide range of Pth- This result 
is similar to a work obtained by E. Bernstein and U. Vazirani J3|. We mention it in Section 9.1. 

Figure [lO] shows a transition about whether quantum computing is available or not for threshold 
probability P t h- Here, let us consider where is a classical searching on the phase diagram of Figure [l0|. 
We assume that we are looking for one item among unsorted 2™ items. If we examine M items from them 
in a classical manner, we can find it with probability P = M/2 n . Now, let us regard M as the number of 
the quantum operations iterated and P as the threshold Pth for the algorithm in classical regime. If we 
give Pth > (and it is not infinitesimal), the classical searching takes M ~ 0(2") and x ^S> x c . Hence, it 
is located far away upward in the non-available region of the quantum algorithm in Figure Uu. 
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xc 




Figure 10: Variation of x c against Pth- A thick solid curve represents perturbation up to the fifth order. 
A thin dashed line shows its tangent at P t h = 1 given by Eq. (|67|). 

On the other hand, if we consider the neighbourhood of P t h = 0, it becomes subtle. The classical 
searching can take small M , and it can approach to the available region of the quantum algorithm. 
Furthermore, in the limit of Pth - » 0, we can expect M t h(p c ) ™ > (but n — > oo) for the quantum 
algorithm, so that behaviour of x c in the neighbourhood of P t h = might be singular. 

From these discussions, we consider that a quantum to classical phase transition of the algorithm is 
described around Pth — in Figure We cannot say anything about it by our approach, because x c 
for P t h — is outside the domain where the perturbation theory is reliable. 

9 Discussions 

In this section, we think about related work obtained by E. Bernstein and U. Vazirani, and how the phase 
error is caused. Then, we give other discussions about our results. 

9.1 Accuracy of quantum gates 

E. Bernstein and U. Vazirani consider accuracy of quantum gates for quantum computation ]^| Jl9|. Let 
us think about a quantum computer which is designed to apply T unitary transformations, Ui, ■ ■ ■ , Ut, 
in succession (T steps) to the initial state \4>o), as follows, 

I^)^|^i)^---^|0t), (68) 

so that \4>t) = Ut\(f>t-i) and (<f)t\(t>t) — (</>o|0o) = 1 f° r t — ,T. On the other hand, we assume that 
it actually applies Ut which is slightly different from Ut to the state because of incomplete accuracy, 

I0o) A H0i}^---^|<M, (69) 

so that \cj)t) = U t \4>t~i), \<Po) = |0o), and {4n\4't) = 1 for t = 1,---,T. (We are considering errors of 
unitary transformations, and it does not cause dissipation to the quantum computer.) 
Defining unnormalised states 

\Bk) = (U t - U t )\<k-i), (70) 

we obtain 

|&) = |0i) + |Pi), 
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\h) = \fo) + \E 2 ) + U2\E 1 ), 



= \4>t) + \E t ) + Ut\E T -i)--- + UtUt-i---U 2 \E 1 ), (71) 

and 

T 

(cfrlfo) >l-J2V(Et\E t ). (72) 
t=l 

Here, we assume that the error of each step is bounded as 

V(Et\E t ) <e fori = 1,---,T. (73) 

We obtain 

\(<f> T \4>T}\ 2 > (1 - Tef = 1 - 2Te + 0(e 2 ). (74) 

Hence, if the error of the unitary transformation at each step is bounded to e, the probability of detecting 
\4>t) that we want as the final state is at least (1 — 2Te). 
On the other hand, from Eq. Ww, we obtain 

(4> t \U T \ct>t-i) = l + {4>t\Et}>l-e, (75) 

and we can rewrite it with density operator as 

tr(pp') > 1 -2e + <3(e 2 ), (76) 

where p — \4>t){<j)t\ and p' = UT\4>t)(4't\Ux- In this paper, we consider the decoherence defined in Eq. (|J). 
Although it is different from the error of unitary transformations in Eq. ((6^), we can obtain 

tv(pp')>l-p, (77) 

and regard it as inaccuracy of operation for e — p/2 at each step. 

If we require |(0t|</>t)| 2 > Pth for a threshold of the probability that the quantum computer gives a 
correct answer, we can obtain 

1 - P th > 1 - |(<M<M| 2 - 2Te, (78) 

as the first order estimation. Substituting e = p/2, and T = 2Mn which is the number of quantum 
gates during the whole process (the number of decoherences caused) into the above, we can obtain 
2Mnp < 1 — P t h- This is similar to the result obtained in Section B, except for a factor. 



9.2 How the phase error occurs 

We give a mechanism which causes the phase error of Eq. (||) for an instance. We can think this error to 
be quite possible for proposed implementations of quantum computation ^ . Let us consider two spin- 
1/2 systems described as 3-component normalised vectors of cr A (qubit) and cr E (environment), whose 
interaction is given by their inner product of k<t a ■ cr E . 

If there is weak external magnetic field along z-direction B — (0, 0, B z ), both of them align themselves 
with z-direction, so that <t a — (0, 0, ±1) and a E = (0, 0, ±1). Hence, we obtain an effective Hamiltonian 
of AH ~ na A af , and a time-evolution operator 



exp(- 


ilAB) 










(00| (01| 


(10| 


(11 




|00) 


fe- ie 










|01) 


e ie 








|10) 





e i0 







111) 


{ 





e -i6 


) 



(79) 
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on the logical basis \ij) = \i A )\i E ) for i,j S {0, 1}, where 9 = (n/Ti)t. 

We assume that the initial state of the systems A and E is given as \f A ) \+ E ), where \<p A ) is an 
arbitrary state of A and \+ E ) = (l/\/2)(|0 B ) + i\l B )). It evolves as follows @, 



= <p A <^ 



where 
and 



pe{+,-} 



M A = (p E \U AE \+ E ), 



\±* 



V2 



(|0 £ )±z|l B ». 



Here, we introduce a 2 x 2 unitary transformation, 



v = = 4= 



)| (i| 

1 10) ( 1 i 

i — i 



and we obtain 



We can rewrite M M as 



E ) = V E \0 E 



VV^ = 7, 



Sli E\ 



V a \l 



M A = 


(o s l 




M A = 




v E ^ 


is given by 






v m u£ E v E = 














" i?t 


" 







R 



u 



AE 



' V 


" 





V 



on the logical basis, where 



Hence, we obtain 



R 



cos a i sin f 
i sin cos 6 



M+ = cos 6I A M_ = -i sin 0cr^ . 
Therefore, the time evolution of the system A is described as 

p A -> p' A = cos 2 6>p A + sin 2 8a z p A a z , 



(80) 

(81) 
(82) 



(83) 
(84) 

(85) 



(86) 



(87) 



(89) 



and it is equivalent to Eq. (||). 

If we assume that each qubit of the quantum computer interacts with an external spin-1/2 particle 
under the weak magnetic field every time interval, our model can give a reasonable description of its 
decoherence. 
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9.3 Other discussions 



From Figure [H], we find x c — 2M t h(x c )np — 0(1) for suitable threshold probability (1/2 < P t h < 1, for 
example). It means that if the error ratio p is smaller than an inverse of the number of quantum gates 
(2Mn) _1 , the algorithm is reliable. If this observation holds good for other quantum algorithms, it can 
serve a strong foundation to realize quantum computation. We cannot investigate a quantum to classical 
phase transition of the algorithm, because it is outside the reliable domain of our perturbation theory. 
For studying it precisely, we may need to construct an exact solvable model of a quantum system with 
decoherence. 
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A Formulas for deriving matrix elements of density operators 

In this section, we collect some formulas that are used for deriving the matrix element of the density 
operator T^ M K 

A.l Formulas of (WR ) k W\0) 

We derive an explicit form of (WRo) k W\0), where |0) is an n-qubit (n > 2) initial state of |0) ® • ■ • ® |0). 
Let us think an n-qubit state of 

|*> =a |0)+ ai Y, \ x )' ( 9 °) 

Using 

W = 7P E (~ 1 T V \V) for Vx G {0, 1}", (91) 

where x ■ y represents an inner product of n-digit binary strings of x, y S {0, 1}™, that is x ■ y — 2~^™=i x iVii 
we can derive 

VUR |*> = W(-a \0)+aiJ2\ x )) 

x^O 

V x^O y 

= * [(- ao + (2™- l)aO|0)-( ao + ai ) ]>>)]. (92) 
Then we introduce a parameter 9 to simplify notations of states fl4|| , 



1 /2 n — 1 

sin6»=^, cos0=W . (93) 

V2" V 2™ 



Using 9, we obtain the following trigonometric formulas, 

:[=Fsin(p+(2 — 1) j = =F sin 9 sin <p + cos 9 cos (p 



/2" V2" - 1 

= cos(ip±9), (94) 

l=[ T cos^-(2"-l)-|=L=] = -sm(v±9), (95) 



V2™ - 1 

^=(±sin<,oH — S = ) = f± cos 9 sin u> + sin 9 cos w] 

fW l V2" - 1 V2 n - 1 
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= ±!^M, ( 96 ) 
V2 n - 1 

1 . since , cos(<p±01 . , 

(±cos^--=^=) - ±W>. (97) 



\/2" V2" - 1 ~ V2" - 1 

From these relations, we can obtain the following results, 

(WiZo)W|0) - cos2l?|0)--==4=5>). (98) 



In general, we obtain 



(WR ) 2k W\0) = (-l) fc [sin(2fc + l)fl|0) + C °^l| ^ E I 1 )] for fc = 0, 1, • • (99) 



(l¥i? ) 2fe+1 ^|0) = (-l) fc [cos2(fc + l)g|0)- Sm ^±i ^El a: )] forfc = 0,l,--, (100) 

A. 2 Formulas of (R o W) k \0) 

From formulas obtained in Appendix [A.l| we obtain 

{R W) 2k \0) = W{WR a ) 2k W\0) 

= {[S in(2, + (2" - 1)^||^]|0) 

+[sin(2fc+ l)^^||±^] E |,) } 

= (-l) fc [cos2fc0|O) + -===E^] forfc = 0,l,---, (101) 
(iioW) 2fe+1 |Q) - (-!)*[- sin(2fc + 1)0)0) + E 1^)] forfc = 0,l,--, (102) 

A. 3 Formulas for summation of trigonometric functions 

When we calculate the matrix element of (0|T^ M ^|0), we often have to sum up trigonometric functions. 
In this paper, we use the following four formulas, which can be proved by the inductive method pjj , 

E(-l)'sin(2* + l)0 = ( ' 1} " Sin2( "t 1)g ten = 0,1,-, (103) 

" COS0 



1=0 



y(-l) z sin2a + l)6» = + sin(2n + 3)9 for n = 0, 1, • • •, (104) 



Z=0 



£(-!)' cob(2Z + 1)« = l + (-l)"cos2(»+l)fl forn = 0il) ... ] (1Q5) 

COS (/ 



Z=0 



n 1 f— 1*1™ 

V(-l)'cos2(Z + l)6» = - + £ — '— cos(2n + 3)9 for n = 0, 1, • • -. (106) 



/ = 
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A.4 Formulas of = (0\(WR o ) l a^(WR o ) k W\Q) 

From Eqs. (|99|), (|100|), (|101[), and (p03), we can derive the following formulas, 



g«(2fc,20 = <0| (V^i^o) 2 '^^ (T^i? ) 2fc V^|0) 

sii 

7^ 



(-l/+ fe [cos2^0| + ^=5:( 



X^)[sin(2 fc + l) g |0> + C0S ^ + 1)g El^ 

= (-l) /+fc [cos 2Z0 sin(2fc + 1)6* - ^ sin 2Z(9 cos(2fc + 1)0] 

for k = 0,1, •••,/ = 1,2,- ■-, (107) 
£ (1) (2fc,2Z + l) = (0\{WR ) 2l+1 a^(WRo) 2k W\0) 

= (-l)' +fe [- sin(2/ + l)0sin(2fc + 1)0 - — - — cos(2Z + l)0cos(2fc + 1)01 

2™ — 1 

for k = 0,1, •••,/ = 0,1,- ••, (108) 
(2k + 1,21) = (0\{WR ) 2l ai l) (WR ) 2k+1 W\0) 

= (-l) i+fe [cos2Z0cos2(fc + 1)0 + — - — sin2Z0sin2(fc + 1)9] 

2™ — 1 

for k = 0,1, •••,/ = 1,2,- ••, (109) 
0^(2*; + 1,21+1) = (0\(WRo) 2l+1 a^(WRo) 2k+1 W\0) 

= (-l) i+fc [- sin(2Z + l)0cos2(fc + 1)9 + — - — cos(2/ + l)0sin2(fc + 1)9] 

2" — 1 

for k = 0,1, •••,/ = 0,1,- ••, (110) 

A. 5 Formulas of (0\(WR o ) k (rP\Vj) and (l/V2)(0|(Wi2 ) V«(l ) - \j)) 

From Eqs. (||), (||) in Section |, and Eqs. (|l0l|) , (|Iol ) in Appendix we can derive the following 
formulas, 

(0\(WR Q ) 2k a^\ Vj ) = (-^[ooB^Ol + ^^XX^^fElf) 
= (-l) fc "!^ fc l = ^-(-l)2- 1 



= -(-l)y F3l sin2fc^ 

(-I) fc sin2fc0„ , . 

-L(0\(WR Q ) 2k aU(\0) - |j)) = -i=(-l)' £ [cos2fc0<O| + -^E ]>>lk«(|0) - |j» 

(-l) fe r , nA sin2fc6> n , . 

= i-^-[cos2fc0-(-l)^- 7 =_], (113) 

-i=(0|(Wi?o) 2fc+ V«(|0)-|J» = -Hg [sin(2fc + 1)0 + (-!)«« + (114) 
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G (h) (l , li, ... , l h ) 



a z (n) CTz (i2) 
x — 



— 



lo 



lh 



Figure 11: Diagram of h errors, G^ h '(lo,h,- ■ ■ ,lh)- 



A. 6 Derivation of the asymptotic forms of matrix elements 

In this section, we derive the asymptotic forms of matrix elements, Eqs. (^4j) and (|5^) introduced in 
Section [?]. 

Comparing Eqs. ( f45f ) and (pl|), we notice that contributions of (0|T^ M \0)/(Mn) h under n — > oo comes 
from only terms in which a z errors occur at different steps and different qubits. The reason is as follows. 
The term of 

| (0| [product of (WR )s and cr 2 s]|0)| 2 (115) 

never exceeds unity, and the number of terms where a z errors occur at the same step or the same qubit 
in (0|T^ M) |0) is at most 0((2M) h n h ~' 1 ). Because we divide (0|T^ M) |0) by (Mn) h , they are eliminated 
under n — > oo, as far as h is finite. 

An asymptotic form of the matrix element with h errors at different steps and qubits, as shown in 
Figure 1 1 , are given by 



lim Q {h) (l ,h,---,lh) 

n — >oo 



lim ({)\{WRo) lh (J z lh) ■ ■ ■ (WRo) h a^ l) (WRo) lo W\0) 



4)£L LW 2 J+ElW sin l 
I cos J c 



((lo + 1)9) 



sin 



(he) ■ 



cos 
sin 



(kO) 



for h = l,2,- 



where lo = 0, 1, • • •, l s = 1, 2, • • • for s = 1, • • • , h, {i s : 1 < i s < n for s = 1, • 
other, [a;J represents the largest integer that is less than or equal to x, 



= Z s (mod 2) e {0,1} for s = 0, • • -,h, 



and 



f(x) for a = 
g(x) for a = l 



(116) 

, h} are different from each 
(117) 

(118) 



We show that this holds for h — 1 in A ppendix A. 4 . We prove Eq. (116) for arbitrary h by the 
inductive method. Let us assume Eq. (116) is satisfied for h. We examine whether it is satisfied for 
(h + 1) or not. 

First, we assume lo — 2m and l\ — 2m (both of them are even). Using Eq. (|3C|), we obtain 
g( h+1 \2m,2m, l 2 ,---,l h +i) 

= (Q\(WRo) lh+1 (T z th+l) ■ ■ ■ (WRo) l2 (J ( z 2) 
~ G^(2(m + m),l 2 ,---,l h+1 ) 



(WRq) 2m af 1 >(WRof m W\Q) 



-V2(-l) m cos(2m + 1)6[-V2 ^ (2k, l 2 , ■ ■ ■ , l h+x ) 



k=0 



+ (0|(^i? )' fc+1 ^ h+l) • • • (WRo) 12 ^^)} 
for m = 0, 1, • • •, m = 1, 2, • • •, as n — > oo, 



(119) 
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where we use lim^^oo cos 9 — 1. (We eliminate one u z operator from the above equation. This technique 
is used in Sections |5| and ^|.) From now on, we use a symbol of '~' as an asymptotic equal sign under 
n —>■ oo for a while. On the other hand, 



(0\(WRo) lh+1 ai ih +^ ■ ■ • (VKi? )M l2) l%) = 0, 



(120) 



which is proved in Appendix |A.7| . Substituting Eq. ( |llq ) for h and Eq. (120) into Eq. (119), and using 
Eq. ( fL03| ), we obtain 

g( h+1 \2m,2m,l 2 ,---,l h+1 ) 

~ sin(2m + 2m + l)ftF(Z 2 , ■ • • , l h+1 ) 

rh — 1 

+2(-l) m cos(2m + 1)0 (-1)* sin ( 2fc + l Wih, ■", k+i) 

( i \rh— 1 Orfi-.fi 

= (-l) m [(-ir sin(2m + 2m + 1)9 + 2cos(2m + l)0^-£ ^-]Hh, • • • , k+i) 

2 COS 9 

~ (-l) m+ ™sin(2m + 1)6> cos 2m9T{l 2r -- ,l h+1 ) asmoo, 



where 



in J 



cos 
sin 



(h+iO). 



(121) 
(122) 



This means Eq. (116) is satisfied for (h + 1) with (Iq, li) = (2m, 2m). 

Next, we assume Iq = 2m and l\ = 2m + 1 (even and odd). From Eqs. (p2|) and (ffit), we obtain 



(WR Q ) 2rh+1 ai ll \WR o ) 2m W\0) 



g( h+1 H2m,2m+l,l 2 ,---,l h+1 ) 

= (0\(WR a ) lh+1 a ( * h + l) ■ ■ ■ (WR ) h a { * 2) 

~ gW((2m + 2m) + l,Z 2 ,---,Z /t +l) 

m — 1 

-\/2(-l) m cos(2m + 1)6[-V2 £ £?W(2fc + 1, Z 2> • • • , l h+1 ) 



k=0 



1 



+ _ (0 |(^i? o) ^ +1(T fe +1 ) . . . (WRo) h 4 2) m - 
for m = 0, 1, • ■ •, m = 1, 2, • • •, as n — > oo. 



(123) 



(We eliminate one cr z operator from the above expression as well as Eq. (119).) On the other hand, 

(0|(W J R o y*+V^+ 1 ) • • • (VFi? )M 42) (|0) - In)) ~ Hh, • • • , h+i) as n - oo, (124) 



which is p rove d in Appendix [A.q . Hence, substituting Eq. ( 116 ) for h and Eq. ( |124| ) into Eq. (123), and 
using Eq.(106), we obtain 

g( h+1 \2 mi 2m+l,l 2l ---,l h+1 ) 
~ {(-i) m+fn cos 2(m + m + 1)0 

m — 1 

+2(-l) m cos(2m + 1)6[J2 (-1)* ™s2(k + 1)9 - -]}T{h, h+i) 

k=0 

= (-l) m {(-l)" l cos2(m + m + l)6» 



1 (-l)m-l 1 

+2 cos(2m + 1)0[- + V o ' cos(2m + 1)9 - -]}T{1 2 , • • • , Z h+1 ) 
1 z cos f z 



(-1) 



m+m+1 



sin(2m + 1)0 sin(2m + 1)&F(Z2 5 • • * , as oo. 



(125) 



This means Eq. flliej ) holds for (h + 1) with (l , h) = (2m, 2m + 1) 
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In the case of (lo,h) = (2m + 1,2m), and (2m + 1,2m + 1), we can show Eq. (116) is satisfied for 
(h + 1) in similar ways. Therefore, we obtain Eq. (116) for h = 1, 2, • • ■ by induction. We pay attention 
that it does not depend on i\, ■ ■ ■ , ih- 

To obtain an asymptotic form of (0|T)[ M - ) |0) /(Mn) h 1 we take intervals of Figure [ll] as lo, • • • , lh-i, 2M— 
Ss=i ' s ' multiply a factor n\/(n— ft.)! to the terms for permutation of a z errors, and sum up them by 
lo, ■ ■ ■ , hi-i, 



lim , , , 

n^oo (Mn) h 



iT (M) |0) 



lim 

n— >oo [Mn) 



^ M-hM-h-lo M-h-{l -\ hlh-2) 

EE-- E 



l =0 h=l Jfc-i=l 



(126) 



Then, substituting Eq. (|116| ), 9 = lim^oo M9, 4> s +i = l s Q for s = 0, • • • , h — 1, 

m-/i e 



lim 

n — >oc ^ — ' 



Zn=0 



and 



lim 

n — >oc 



Z-<i=o r®-> <frt 

E ' ' 



ls=l 



for s = 1, • • • , h — 1, 



(127) 



(128) 



into Eq. fll26j ), we obtain Eqs. (|5J and (|56j). 

A. 7 Formulas of (0\(WR o ) lh cr^ ■ ■ ■ (WR ) h a^\r] io ) 
In this section, we show 

(0\(WR o ) lh a^ ■ ■ ■ (WR ) h ai^\r, l0 ) = 0, (129) 
for ft = 1, 2, • • •, where Z s = 1, 2, • • • for s = 1, • • • , h, an d \i s : 1 < i s < n for s = 0, 1, • • • , ft,} are different 



from each other. For ft, = 1, it is proved in Appendix A. 5 
To prove Eq. (129) for arbitrary ft, we need to show 

{OKWRoY^i^ ■ ■ ■ (WRo) h a^\r, ju ..., js ) = 

for • • • , ih, ji, ■ ■ ■ , j s } that are different from each other, 
(0\(WR o ) lh <ji^ ■ ■ ■ (WR ) h ai^(\j ,ju---iJs) - \n~~Ts)) = 

for {*!,••■, ih,jo,ji, ■ ■ ■ , js} that are different from each other, 



where 1 < j <«,•••, 1 < j s < n, 



\Vh, 



V2™ 7 



E i*>. 



and 



.y'o-.y i. ••••./- i ')•••() l ()•••() l ()•••<) . 

T T T 

JO jl js 



(130) 
(131) 

(132) 
(133) 



From similar calculations in Appendix A. 5, we can show Eqs. (130) and (131) are satisfied for ft = 1. We 
prove them for arbitrary ft by the inductive method. 
Let us derive an explicit form of 



(WR ) , <T^\f th ,... J .) for « = 1,2,---, 



(134) 



27 



which appears in Eq. (130). At first, we consider the case of s = 1. We can obtain 



yj = 1 yj = 1 

Vi=° Vi =1 



= \Vj) - ^2\Vj,i) fori^j. 
Then, we derive an explicit form of (WRo) \T]i,j)- In the case of I = 1, we obtain 

1 



WR o \r h , j ) = -(\0)-\i)-\j) + \i,j)). 



(135) 



(136) 



For Z = 2, we obtain 



(WR ) 2 \n 



1,3 I 



: __ W (| ) + |i)+| i )-|i,j)) 

= -W\Q) + \thj). 

Hence, we obtain 

l-i 

{WR Q f l \m,j) = -Y,{WR Q f k W\0) + \ m , j ) fori = 0,l,---, 

fe=0 

l-i 

Y,(WRo) 2k+1 W\0) + -(|0) - \i) - \j) + for I = 0, 1, 



(WR ) 2l+1 \mj) 



(137) 

(138) 
(139) 



k=0 



1-1 



From Eqs. (g8|), ( |135| ), and ( |138| ), we obtain 

(WRo) 2l a^\ Vj ) = (WRo) 2 %)-V2(WR ) 2 %J 

i-i 

= -V2j2(WRo) 2k W\0) + \ VJ ) + v / 2^(Wi? ) 2 'W|0) - V%-,i> 

= \Vj) ~ V2\Vj,i) 

= <T®\r/j) for l = 0,l,---. 

Furthermore, using Eqs. (^9|) and (|139| ), we obtain 

(W^of+VMh,) = WRo<ri%j) 

= WR {\ m )-V2\ m ^) 

= for/ = 0,l,--. 

In general, from similar calculations above, we can obtain 

(WRo) 2l a^\ Vju ... Js ) = \r) ju :, jB ) - ^n^J &r / = 0, 1, • • •, 
{WR ) 2l+1 a^\r lju ...^) 

= i E {-\r + - +a -\i, ai j u ...,asjs) 

* ae{0,l}» 

= -j= (-l) Q2+ '" +Q3 (N,0,a 2 j2, • • -,a s j s ) - \i,ji,a 2 j 2 ,-- -,a s j s )) 



(a 3) — ,a s )E{0,l}" 

for 1 = 0, 



(140) 



(141) 



(142) 
(143) 



(144) 
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where 



|i,aiji,-",a fl j s ) =|0-0 10-0oi0-0a, 0-0). 

T T T 

i h js 



(145) 



Here, we notice that, if we relabel in dexe s, the right side of Eq. ( 144 ) is written as a sum of (\jo, ji, • • • , jV 
lilj ' ' ' )Jr))j which appears in Eq. ( |l3l| ). 
Then, we derive an explicit form of 



(W^ )VW(|j ,ii, • • • , Js) - • • • , j s )) for J = 1,2,- • •, 



(146) 



included by Eq. (131). Here, we can assume jt = t for t = l,---,s, and s < i < n, without losing 
generality. Because {i,jo, • • ■ > js} are different from each other, ai^ does not have an effect on the state. 
Hence, we obtain 



(WRo) 1 <t®(\ 1,- ••, 1,(^0) - |0,l } ---,l,0,---,0» 

a+l n— s— 1 s n— s— 1 

V--,i,o, — .oMO^r-^tV^o) 

n — s — 1 
Xl,- 



V2^ 



S+l 71 — S — 1 S 

' (-l) 1+Sl+ "' +a: il^i,---^n-i) for i= 1,3, 



for I = 0, 2, • • • (even) 
(odd) 



(147) 



cG{0,l}" 



From this result, we find that Eq. (146) is described as a sum of (|jo> ji, • ■ • , js) — |ji> ■ ■ 1 , js)), which 
appears in Eq. (131). 



Here, we assume Eqs. ( |130D and (131) hold for h. From Eqs. (U42), (144), and (147), we can resolve 
Eqs. (|130|) and (131) for (h + 1) to terms that contain h errors of cr z , so that we can prove them to be true 



for (h + 1). Hence, they are proved for any h by the inductive method. Therefore, we obtain Eq. Ql29| 
for arbitrary h. 

A.8 Formulas of (0\(WR ) lh a^ ■ ■ ■ (WR ) h a^\\0) - 
In this section, we show 

(0\(WR o ) lh a^ ■ ■ ■ (WR ) h 4 h) (\Q) - |*o» ~ Hh, •••,//.) as n -> oc, 



(148) 



for h = 1, 2, • • •, where Z s = 1, 2, • • • for s — 1, • • • , h, {i s ■ 1 < i s < n for s = 0, 1, • • • , h} are different 
from each other, and the function T(l\, ■ ■ ■ , lh) is defined in Eq. (122) of Appen dix [A.6| . Eq. ( |148| ) is used 
for the inductive method in Appendix |A.6| . It is shown for h — 1 in Appendix A. 5. We prove it for any 
h by the inductive method. 



We assume Eq. (148) is satisfied for (h — 1). Let us consider the following equation, 

(W J R )V«(|0> - |J» = -(WRoY-'WiW + m 

= -(WR ) l ~ 1 (2W\0) far Z = 1,2,- ■■ and i^j. 



(149) 



Here, we assume £ = 2m + 1 (odd). From Eq. ( |28[) , we obtain 

(H/i?o) 2m+1 CT«(|0) - |7)) = -2{WR Q ) 2m W\0) + V2(WRo) 2m \Vj) 

m— 1 

= -2(WR Q ) 2m W\0) -2^2 (WR ) 2k W\Q) + V2|jft) 



-2 ^(Wi?o) 2fc W|0) +72^). 



fc=0 



Using Eq. (129), we obtain 

rn 

(0\(WR o ) lh ai^ ■ ■ ■ (WR ) h -ai^(WRo) 2m+1 a^\\0) - = -2^ g^(2k, Z 2 , • • • , k) 



(150) 



(151) 



fc=0 
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Then, we consider the following fact. We use Eq. (148) for obtaining the asymptotic form of Eq. ( |116[ ) 
by induction in Appendix A. 6. Hence, we can assume that Eq. ( |116| ) is true for (h — 1). Therefore, we 
can require Eq. (151) as 

(OKWRoY^i 1 ^ ■ ■ ■ {WRo) l2 (ri i2) (WRo) 2m+1 a^dO) - 

rn 

~ -2^(-l) fe sin(2fc + l)9F(l 2 , ■■-,l h ) 

k=0 
\m+l , 



(-l) m+1 sin 2(m + l)0.F(Z 2 ,-- ■ ,l h ) 



T{2i 



1,/a 



as n 



oo, 



(152) 



where we use Eq. (103) and sin2(m+ 1)6* ~ sin(2m+ 1)9. 

Hence, Eq. (148) is satisfied for h in the case that I is odd. For I = 2m (even), we can give a similar 
discussion and prove it. Therefore, we obtain Eq. (|148|) for arbitrary h. 



B Notes for numerical calculations 

In this section, we take some notes about numerical calculations of higher order perturbations. 

First, we calculate the asymptotic forms of the forth and fifth corrections for density operator. Using 
the rules of Eqs. (p3) and (pq), we obtain 



r.,™ 1 15- 169 2 in 5(3 + 326 2 ) . ^ 

F ^ Q) = 48 + T2288e^ C ° s4e ' 4915263 Sm 49 

= le 2 + o(e 3 ), 

^ 1 45 + 7206 2 - 2566 4 , n 3 + 326 2 + 2566 4 . ^ 

F 5 (Q) = 1 -. cos 46 = sin 46 

K 1 240 19660806 4 5242886 s 

= I ^6 2 + 0(6 3 ), (153) 

where i*\(0) is defined in Eq. (|6l|). From the definition of Eq. (|60|), we can obtain 64(6) and £5(6), the 
forth and fifth coefficients of Eq. ( |59| ) . 

Next, we consider the region of x where the perturbation theory up to the fifth order is valid. To 
investigate it, we examine the sixth order correction, which is given by 

1 315+ 16806 2 - 2566 4 ^ 7(15 + 2566 4 ) . ^ 

F 6 (Q) = 1 -. cos46 ^ sin 49 

v ' 1440 235929606 4 314572806 5 

= ^ + °( 03 )- ( 154 ) 
From numerical calculations, we obtain 

< ^C 6 (6) < 1.62 x 10" 4 , (155) 
6! 

where < 6 < (7r/4). Hence, if we limit x to 

0<x< 1.35, (156) 

it is bounded to 

< ^C 6 (9)i 6 re 9.79 x 10~ 4 < 10" 3 . (157) 
6! 
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